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Abstract 

The  paper  presents  a  set  of  two-dimensional  molten  carbonate  fuel  cell  mathematical  models,  considering  the  electrical,  mass  and  heat  transfer 
characteristics  of  a  molten  carbonate  fuel  cell.  A  dynamic  simulation  model  was  built  under  the  VC++  environment.  The  steady  state  and  dynamic 
simulation  results  for  a  cross-flow  molten  carbonate  fuel  cell  were  obtained  using  the  simulation  model.  The  results  are  appropriate  for  the  design 
and  operation  of  molten  carbonate  fuel  cells. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

A  fuel  cell  is  considered  as  a  good  candidate  for  the  next 
generation  power  source.  Fuel  cell-based  power  plants  convert 
chemical  energy  directly  to  electricity  without  the  requirement 
for  conversion  of  energy  into  heat.  In  addition,  NO*  and  SO* 
emissions  are  greatly  reduced  in  comparison  with  fossil  fuel- 
based  generation  systems  [1,2]. 

Among  the  various  fuel  cell  types,  the  molten  carbonate  fuel 
cell  (MCFC)  is  a  very  promising  type,  which  is  now  about  to 
become  commercially  used.  MCFCs  belong  to  the  high  tem¬ 
perature  fuel  cell  class:  the  operating  temperature  range  is  about 
600-700  °C,  in  which  the  electrolyte  has  a  good  ionic  conductiv¬ 
ity  [3].  The  MCFC  is  considered  to  be  a  second-generation  fuel 
cell.  It  has  high  efficiency  and  low  emissions  of  NO*  and  SO*. 
The  MCFC  has  been  proposed  in  recent  years  as  an  extremely 
efficient  solution  in  power  generation. 

The  distribution  parameters,  such  as  the  reaction  temperature, 
current  density  and  gas  species  concentration,  have  important 
effects  on  the  operation  of  the  fuel  cell.  So  they  will  have  a  very 
important  significance  in  predicting  the  distribution  of  temper¬ 
ature,  current  density  and  species  concentration  in  a  fuel  cell 
stack.  Considering  the  electrical  and  heat  characteristics  of  a 
fuel  cell,  the  paper  presents  a  two-dimensional  dynamic  mathe¬ 
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matical  model  of  the  MCFC.  The  steady  and  dynamic  simulation 
results  of  a  cross-flow  molten  carbonate  fuel  cell  were  obtained 
using  the  simulation  model. 

2.  Working  principles  of  a  MCFC 

The  working  principle  of  a  fuel  cell  is  the  simultaneous 
reduction  of  an  oxidant  and  oxidant  of  a  fuel.  Any  oxidant 
gas  containing  oxygen  and  carbon  dioxide  is  supplied  to  the 
cathode.  There  the  oxygen  reacts  with  carbon  dioxide  forming 
carbonate  ions.  These  ions  are  transported  towards  the  anode 
where  they  are  involved  in  the  oxidation  reaction  of  hydrogen. 
At  the  anode,  water  and  carbon  dioxide  are  formed.  The  anode 
reaction  produces  electrons,  while  at  the  cathode,  electrons  are 
consumed.  The  electrons  pass  from  the  fuel  electrode  to  the  oxi¬ 
dant  electrode  via  an  external  circuit,  and  the  electrical  circuit  is 
completed  by  ions  crossing  through  an  electrolyte.  The  layout 
of  the  MCFC  is  showed  in  Fig.  1. 

At  the  fuel  cell  electrodes,  the  following  electrochemical 
reactions  occur:  [5] 

At  the  cathode: 

j02  +  C02  +  2e“  =►  C03-2  (1) 

At  the  anode: 

H2  +  C03~2  =►  H20  +  C02  +  2e“  (2) 
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Nomenclature 

E 

cell  potential  (V) 

F 

faraday  constant  (c  mol-1) 

h 

enthalpy  (kJkg-1) 

i 

local  current  density  (Am-2) 

K 

equilibrium  constant 

h 

stoichiometric  reaction  number 

An 

molar  flow  rate  change  (molm-3  s-1) 

A  P 

pressure  drop  (bar) 

a 

Generation  heat  (kJ  m-3  s-1) 

Q 

reaction  heat  (kJ  kg-1) 

^local 

internal  resistance  (Q  m2) 

R 

mass  rate 

AS 

entropy  change  of  the  reaction  (J  mol-1  K-1) 

T 

temperature  (K) 

V 

working  voltage  of  fuel  cell  (V) 

AV 

local  voltage  losses  (V) 

W 

power  (kW) 

y 

mole  concentration  (mol  mol-1) 

0 

polarization  resistance  (£2  m2) 

Subscript 

A 

anode 

C 

cathode 

E 

electrode 

MCFC 

molten  carbonate  fuel  cell 

SI 

up  separator 

S2 

down  separator 

sh 

shift  reaction 

heat.  Meanwhile  there  is  mass-transfer  between  the  gas  of  anode 
or  cathode  and  the  electrolyte,  which  forms  a  mass  flow  along 
with  a  heat  flow.  And  when  the  current  flows  through  fuel  cell, 
this  produces  resistance  heating  due  to  the  various  resistances. 
The  heat  transfers  between  the  gas  and  solid  parts  by  conduction, 
convection  and  radiation.  It  will  cause  a  change  in  the  fuel  cell’s 
temperature  field.  At  the  same  time,  the  temperature  distribution 
also  affects  the  electrochemical  reaction. 

It  can  be  seen  that  the  thermal  physical  and  electrochemical 
processes  are  interactive  and  interdependent  in  the  fuel  cell.  The 
distribution  of  temperature  has  a  great  effect  on  the  electrochem¬ 
ical  reaction.  So  the  analysis  of  heat  generation  and  transfer  in 
a  fuel  cell  is  important. 

3.  Mathematical  model  for  the  MCFC 

Several  assumptions  as  following  are  applied  in  this  paper: 

(1)  The  gas  flow  is  in  smooth- walled  channels  of  constant  cross 
section. 

(2)  Cell  voltage  has  no  distribution  and  is  constant. 

(3)  The  electrolyte,  anode  and  cathode  are  considered  as  a 
whole,  they  have  the  same  temperature  and  character  at  the 
same  position. 

(4)  The  fuel  cell  wall  is  completely  insulated,  and  heat  loss  is 
neglected. 

(5)  The  cell  is  with  cross-flow  in  which  fuel  gas  is  in  the  x- 
direction  and  oxidant  gas  is  in  y-direction. 

(6)  The  calculated  infinitesimal  is  orthogonal. 

x  =  1 ,  . . . ,  j,  . . . ,  N,  y  =  1 ,  . . . ,  k,  . . . ,  M 


The  overall  electrochemical  reaction: 

H2  +  ^C>2  (g)  +  C02(c)(g)  =  H20(g)  +  C02(A)  (3) 

There  are  many  processes  in  the  fuel  cell,  including  heat 
generation  and  transfer,  mass-transfer,  oxidation,  ionization,  ref¬ 
ormation  and  so  on.  Electro-chemical  reactions  can  produce 


3.1.  Power- generation  characteristics 

The  electrochemical  reactions  are  considered  to  be  very  fast, 
and  modeled  in  a  steady-state  form. 

The  cell-unit  performance  of  the  electricity  generation  is 
expressed  as:  [6] 


^  2e  Up  seperator 

H?  1  T  ca  T  H'° 

'  anode 

A 

T 

CO:~  electrolyte 

t 

^  cathode 

Ji  '0 

CO,  -O, 

2e  — ►  Down  seperator 

V  —  E  A  Vohmic  AVanode  AVcathod  (4) 


where  V  (V)  is  the  cell  operating  voltage,  E  (V)  the  reversible 
local  cell  potential,  AV0hmic  the  local  voltage  losses  due  to  the 
internal  ohmic  resistance,  AVanode  and  AVcathod  are  the  local 
voltage  losses  due  to  an  overpotential  at  the  anode  and  cath¬ 
ode,  respectively,  R0hm  the  internal  resistance,  and  77anode  and 
^cathode  are  the  polarization  resistance  at  the  anode  and  cathode, 
respectively. 

According  to  the  Nemst  equation: 


E=  E°  + 


RT  ]  yn2  x  }’of  x  yco2(c) 

- In - - - 

2E  yH2o  x  37co2(A) 


RT 

H - In  P 

IF 


where  E°  is  the  standard  potential,  P  (bar)  is  the  gas  pressure. 


2 

^local  =  (^ohm  T  0 anode  T  ^cathode)  m  ) 


Fig.  1.  MCFC  schematic. 
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V  —  E  i  x  (^ohm  T  0 anode  T  ^cathode)  (V) 


(6) 


Internal  cell  resistance  is  usually  expressed  by  means  of  an 
Arrhenius  equation  as  a  function  of  temperature  [4] : 


Rohm  =  0.5  x  10  4  x 


exp 


,  1  1 

3016 - 

T  923 


T 


co2  "H20 


^cathode  =  7.505  X  10  10  X  exp 


(  9289  \  n_o.43  n-o. 


\  T  )  <  h- 


09 


CO2 


(7) 


rianode  =  2.27  x  1(T9  x  exp  (  ^  )  P^A2  Prop  P»}n  (8) 


(9) 


It  can  be  seen  that  the  distribution  of  temperature  has  direct 
effect  on  local  potential  and  resistance  from  above  formulas. 


3.2.  Conservation  of  chemical  species 

During  electrochemical  reactions,  two  mole  electrons  are 
generated  with  one  mole  H2  at  the  anode,  and  consuming  a  half 
mole  O2  and  one  mole  CO2  at  the  cathode. 

The  consumption  or  liberation  flow  rate  in  (molm-2s-1)  of 
the  gas  species  H2,  N2,  CO,  CO2,  O2  and  H2O  is  related  to  the 
local  current  i  (A  m-2)  generation: 

i 

m  =  n  •  —  (10) 

IF 

m  =  (wh2,  ^n2,^co,  ^co2,  ^o2,  ^h2o)(1) 

where  F  is  the  Faraday  constant  (96493  c  mol-1);  i  the  local  cur¬ 
rent  density,  and  h  stand  for  the  stoichiometric  reaction  number 
of  individual  gas  species. 

In  the  fuel  gas,  the  water-shift  reaction  takes  place: 

H2  +  C02^  H20  +  CO  (11) 

The  composition  change  Arshin  (molm-3  s-1)  due  to  the  shift 
reaction  is  determined  by  the  local  fuel  gas  concentration.  The 
reaction  is  double  direction. 

These  chemical  processes  are  very  fast  on  a  time  scale  of 
10  jxs.  They  are  fast  processes  compared  to  the  heat  flow  in  a  fuel 
cell.  Moreover,  volume  effect  of  infinitesimal  was  considered 
very  small.  Therefore,  it  is  reasonable  that  the  electrochemical 
processes  are  modeled  in  steady-state  form. 


3.2.2.  Conservation  of  chemical  species  for  the  cathode  gas 
For  the  cathode  gas 

(Fk- 1  +  ^2  x  yk  =  Fk- 1  x  yk~  1  +  Ri  (13) 


where: 

E*< 


3/ 

—  x  AT  x  AX 
4  F 


Ri  = 


IF 


x  ^cathode  x  ^surface 


^cathode (N2,  O2,  CO2,  H2O)  =  (0,  ,  -1,  0) 


3.3.  Conservation  of  energy  for  fuel  cell 

The  study  of  heat  generation  and  transfer  in  fuel  cell  is  most 
important  because  it  has  a  close  relation  to  the  electrochemical 
reaction  rate  and  power-generation  characteristics. 

The  heat  generation  in  the  cell  unit  arises  from  three  aspects: 
the  electrochemical  reaction,  the  cell  losses  and  the  shift  reac¬ 
tion. 

The  local  heat  generation  in  the  cell  unit  arises  from  the  elec¬ 
trochemical  reaction  and  cell  losses  qQ  (kJ  m-3  s-1): 

i  x  (— (T  AS/2F)  +  ((£  -  V)/1000)) 

<7e  =  - : -  (14) 

Cell 

where  AS  is  the  entropy  change  of  the  reaction,  tce\\  (m)  is  the 
thickness  of  the  cell  unit. 

The  heat  generation  due  to  water-shift  reaction  is  qs: 

qs  =  A //Shift  x  A«co  (kJ  m-3  s-1)  (15) 

where,  A//Shift  is  the  enthalpy  change  of  the  shift  reaction  and 
Anco  (molm-3  s-1)  is  the  molar  flow  rate  change  in  the  fuel 
gas  volume  unit. 


3.3.1.  The  energy  equation  for  fuel  and  oxidant  gas 
The  general  equation  is  expressed  as: 


d(pxh)  f  -  kVh\ 

+  div  (  p  x  V  x  h - — —  =  Oh 


dr 


C p  ) 


(16) 


The  up-wind  difference  equation  of  the  fuel  gas  in  X  direc¬ 
tion: 


3.2.1.  Conservation  equation  of  chemical  species  for  the 
anode  gas 

For  the  anode  gas 

(Fj~  i  +  X, R,)  x  yj =  Fi- 1  x  v-i  +  R<  v'F) 

where: 


AX  x  AT  x  Ha  x 


P  dT 

x  Cpj  x 


A 


RxT 


1 


dr 


=  Fj-\  x  Cpj- 1  x  Tf_x  -  (F j-\  +  R‘)  x  CPj 


xTf  +  Aea  x  KEa  x  (Tf  -  Tf)  +  ASa  x  Ks a 


-  ->  l 

Ri  —  An  shift  X  ^shift  X  V  +  — —  X  nanode  X  Asu 

2F 


^ shift (h2,  CO,  C02,  H20)  =  (1,  -1,  1,  -1) 


x(r/sl  -  Tf)  +  qs  x  AX  x  AF  x  HA  +  ^  x  A 

xCp  x  Tf  -| - x  A  x  Cp  x  Tf 

y  J  2F  j 


^ anode (H2,  CO,  C02,  H20)  =  (-1,  0,  1,  1) 


i 

IF 


x  Ax  Cpu2  x  Tf 


(17) 
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The  up-wind  difference  equation  of  the  fuel  gas  in  T  direc¬ 
tion: 

P  d  TP 

AX  x  AY  x  Hcx  -  x  Cpk  x  — ^ 

RxTk  y  dr 

=  Fk~ i  x  Cpk- 1  x  7f_j  -  (Fk_ i  +  ^  /?,•) 
xCpk  xr[  +  Aec  x  ^ec  x  (,7/C  —  7/C)  +  Asc 

xASC  x  (7f2  -  if)  -J-xAXxAYxCp 

4  F 


xf;L  - 


x  AXx  AY  xCpxTl 


(18) 


The  difference  equation  for  down  separator: 


AX  x  AT  x  //s  x  ps  x  CpS  x 


=  AX  x  AT  x  Hs  x  As  x 


dTS2 


7^S2  _  9  V  7~'S2  t-'S2 

9+U  z  x  El  +  Lj~  C 


rrS2 

ljM  1 


9  v  y-S2  i  j-S2 

z  x  fjA  +  E*~ 1 
AF2 


AX2 


-  AA  x  AF  x  XSC 


x(^ 


E*)  + 


CTB  X  AX  x  AF  x  (7f?  - 
(1Aec)  +  (1As2)- 


tS24a 


(22) 


J.X2.  77ze  energy  equation  for  fuel  cell  (including: 
electrolyte,  anode  and  cathode )  and  separators  (including: 
up  and  down) 


AX  is  length  of  infinitesimal,  AT  is  width  of  infinitesimal,  He, 
Ha,  Hq,  Hsi,  Hs2  are  the  thickness  of  electrolyte  and  electrode, 
the  height  of  anode,  the  height  of  cathode,  the  thickness  of  up 
separator  and  the  thickness  of  down  separator,  respectively. 


dT  (  d2T  d2T  d2T\ 

"xC-x  37-(,A'x  a?  +  A’x  3F+a«x  3?J=« 

(19) 


The  difference  format  of  two-dimensional  diffusion  item  is 
FTCS. 

The  difference  equation  for  fuel  cell  (electrolyte,  anode  and 
cathode): 


AX  x  AT  x  He  x  pe  x  Cpe  x 


dTH 


4.  Simulation  and  analyses 

The  fields  of  temperature,  current  density  and  species  con¬ 
centration  interfere  with  each  other.  To  calculate  the  steady  state 
condition,  the  fields  of  temperature,  current  density  and  species 
concentration  should  be  iterated  until  the  results  converge.  We 
can  obtain  the  distribution  of  temperature,  current  density  and 
gas  species  concentration  of  fuel  cell  at  the  steady- state  point. 
When  the  boundary  conditions  change,  the  fields  of  tempera¬ 
ture,  current  density  and  species  concentration  change  with  time 


=  AX  x  AT  x  He  x  Ae  x 


Tf+l  k-2  xTfk  +  7*  k  Tfk+l  -2  xTfk  +  T* 

J,K  J  1,/C  .  JO CT1  JO C  J,R  1 


AX2 


AF2 


+  qQ  x  AX  x  AT  x  He 


-  AX  x  AT  x  XEA  x  (Tfk  -  T$k)  -  AX  x  AT  x  XEC  x  (Tfk  -  7^)  +  AX  x  AT  x  —  x  Cpc02  x  T*fk  +  AX 


xAT  x  —  x  CD09  x  Tfk  +  AX  x  AT  x  —  x  Cdh9  x  Tyk  —  AX  x  AT  x  — 

4 p  pU2  J*k  2 F  p  2  J,k  2 F 


<rB  X  A  A  x  AF  x  (Tf4  -  Tf4) 
(l/fiEA)  +  (l/£sl)  —  1 


xCpc02  X  Tjk  -  AX  x  AY  x  —  x  CpCH2o 

crB  x  AX  x  AF  x  (T%  -  Tf24) 

(l/eEc)  +  (l/£s2)-l 


rj^E 

X  TJ.k  ~ 


(20) 


The  difference  equation  for  up  separator: 


AX  x  AT  x  Hs  x  ps  x  Cps  x 


dTsl 

aij,k 


=  AX  x  AT  x  Hs  x  As  x 


7^S1  _  9  v  T"S1  _i_  7^S1 

1  j+l,k  1  j,k  '  1  j-l,k 


7^S1 
lfk+ 1 


A  A2 


2  x  +  rg_i 
AF2 


aB  x  A  A  x  AF  x  (Tf4  -  Tf4) 
(1/sea) +  (l/£si)-  1 


x  KSA  x  (Tf  -  Tf)  + 


-  AX  x  AY 


(21) 
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Fig.  2.  Dynamic  response  of  the  current  density. 


correspondingly.  They  can  then  reach  another  steady- state  after 
a  period  of  time. 

4.1.  Dynamic  results  and  discussions 

In  our  case  study,  the  typical  operating  conditions  were  as  fol¬ 
lows:  fuel  gas  flow  rate  and  inlet  temperature  0.017746  mol  s-1 
and  873  K;  supplied  fuel  is  composed  of  10%  CO,  68%  H2, 12% 
CO2,  10%  H2O  (mole  fraction);  oxidant  (air)  flow  rate  and  inlet 
temperature  are  0. 1412  mol  s-1  and  873  K,  the  oxidant  gas  is 
composed  of  14.7%  02,  54.6%  N2,  30%  C02,  0.7%  H20  (mole 
fraction),  respectively.  The  temperature  of  the  components  in 
the  fuel  cell  is  assumed  to  be  874  K. 

The  reaction  gas  flow  into  the  fuel  cell  at  time  of  0  s.  It  reaches 
steady-state  condition  at  a  time  of  100  s.  At  a  time  of  250  s, 
the  flow  rate  of  the  fuel  gas  and  the  oxidant  gas  step  increase 
20%  at  the  same  time.  After  about  400  s,  fuel  cell  achieves 
another  steady  state.  We  can  obtain  the  main  parameters’  two- 
dimensional  field  of  fuel  cell. 

Five  representative  points  including  [7  =  5]  [7  =  5],  [7=10] 
[7=  1],  [/=  1]  [7=  10],  [/=  10]  [7=  10],  [/=  1]  [7=  1]  were  selected 
as  examples.  The  dynamic  character  of  fuel  cell  can  be  seen  from 
the  dynamic  response  of  the  individual  parameters  at  five  points. 

Fig.  2  shows  the  dynamic  response  of  the  current  density  of 
each  point.  In  the  initial  stage,  the  current  density  is  increasing 
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Fig.  3.  Dynamic  response  of  the  temperature. 


because  the  electrochemical  reaction  is  strengthening  gradually. 
The  average  current  density  of  the  fuel  cell  is  1440  A  m-2  when 
it  reaches  the  steady  state  after  100  s.  At  250  s,  the  flow  rate  of 
the  fuel  gas  and  oxidant  gas  stepwise  increase  by  20%.  It  can 
be  seen  that  the  current  density  then  increases  quickly  due  to 
the  increase  of  species  concentration  and  then  decreases  slowly 
under  the  influence  of  the  temperature.  The  current  density  incre¬ 
ment  at  each  point  is  different.  The  average  current  density  of 
the  fuel  cell  is  15 10 Am-2  when  it  reaches  the  new  steady 
state.  The  output  power  of  the  fuel  cell  is  increased  by  about 
5%. 

Fig.  3  shows  the  dynamic  response  of  the  temperature. 
From  this  figure,  it  can  be  seen  that  the  temperature  of  the 
fuel  cell  is  increasing  in  the  initial  stage  because  of  the 
heat  generation  of  the  electrochemical  reactions  and  the  resis¬ 
tance.  The  increase  of  temperature  accelerates  the  electrochem¬ 
ical  reactions.  So  the  current  density  is  increasing  and  the 
species  concentration  is  decreasing.  The  fields  of  tempera¬ 
ture,  current  density  and  species  concentration  reach  a  steady 
state  after  100  s.  But  the  temperature  of  each  point  is  differ¬ 
ent. 

After  250  s,  the  flow  rate  of  the  fuel  gas  and  oxidant  gas  step¬ 
wise  increase  by  20%.  The  average  current  density  then  increases 
by  about  7%.  The  electrochemical  reaction  is  accelerated.  The 
current  density  increment  of  the  different  points  is  different  and 
the  heat-exchange  is  uneven.  So  some  points’  temperature  is 
increased,  others  decreased.  After  400  s,  the  field  temperature 
reaches  another  steady  state. 
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Fig.  4.  Two-dimensional  field  of  current  density. 


Fig.  6.  Two-dimensional  H2  concentration  field. 


4.2.  Steady -state  results  and  discussion 

The  two-dimensional  distribution  of  current  density  is  illus¬ 
trated  in  Fig.  4.  It  can  be  seen  that  the  point  of  maximum  current 
density  is  located  near  the  inlet  of  the  fuel  gas  flow  and  the 
outlet  of  the  oxidant  gas  flow.  The  electrochemical  reaction  in 
this  region  is  most  intense  because  of  the  highest  H2  concen¬ 
tration  of  the  fuel  gas  and  highest  temperature  of  the  oxidant 
gas.  The  trend  in  current  density  is  decreasing  in  the  direction  of 
the  fuel  gas  flow  and  increasing  in  the  direction  of  the  oxidant 
gas  flow. 

Fig.  5  reflects  the  two-dimensional  temperature  characteris¬ 
tics  of  the  electrolyte.  The  area  of  the  inlet  of  the  fuel  gas  flow 
and  the  outlet  of  the  oxidant  gas  flow  has  the  highest  temper¬ 
ature  because  of  the  most  intense  electrochemical  reactions.  It 
can  also  be  seen  that  the  slope  of  temperature  is  sharp  in  the 
direction  of  the  oxidant  gas  flow.  The  flow  rate  of  the  oxidant 
gas  is  much  greater  than  that  of  the  fuel  gas;  the  oxidant  gas  has 
a  very  strong  cooling  effect  for  the  electrolyte.  Therefore,  the 
oxidant  gas  flow  rate  can  be  used  to  control  the  temperature  of 
the  fuel  cell. 

Fig.  6  shows  the  two-dimensional  H2  concentration  field 
of  the  fuel  gas.  From  Fig.  6,  it  can  be  seen  that  the  con¬ 
centration  of  H2  decreases  in  the  direction  of  the  fuel  gas 
flow  due  to  the  consumption  by  the  electrochemical  reaction. 


10  * 

Fig.  5.  Two-dimensional  temperature  field  of  the  electrolyte. 
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Fig.  7.  Two-dimensional  CO2  concentration  field  of  the  oxidant  gas. 


Moreover,  the  temperature  of  the  oxidant  gas  increases  in  the 
direction  of  the  oxidant  gas  flow,  so  that  the  H2  concentration 
of  the  fuel  gas  decreases  in  the  direction  of  the  oxidant  gas 
flow. 

The  CO2  concentration  of  the  oxidant  gas  decreases  in  the 
direction  of  the  oxidant  gas  flow  because  of  the  consumption 
by  the  electrochemical  reaction,  which  can  be  seen  in  Fig.  7. 
In  addition,  the  area  of  the  inlet  of  fuel  gas  flow  and  output  of 
oxidant  gas  flow  has  the  greatest  current  density,  so  the  concen¬ 
tration  of  CO2  in  this  region  is  lower. 


5.  Conclusions 

A  two-dimensional  dynamic  mathematical  model  of  the 
MCFC  in  this  paper  can  reflect  the  distribution  of  temperature, 
current  density  and  species  concentrations  effectively.  More¬ 
over,  we  can  obtain  the  steady  state  and  dynamic  characteristics 
of  a  cross-flow  molten  carbonate  fuel  cell  using  our  simulation 
model. 

It  can  be  seen  that  the  distribution  parameters  such  as  tem¬ 
perature,  current  density  and  species  concentration  are  affected 
interactively.  This  has  a  very  important  influence  on  the  safety 
and  efficient  operation  of  the  fuel  cell. 
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